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Abstract 

The conditions for the HBT interferometry from stars are compared to those en- 
countered in heavy ion reactions. As a consequence the formalism using wavepackets 
is developed. The modifications relative to the standard formalism with plane waves 
are presented, and the effects of the residual Coulomb interactions between charged 
pions and with the charged source are calculated. 

1 Introduction 

Originally the interferometric method of Hanburry Brown, Twiss (HBT) was devel- 
oped in astronomy and used to determine the angular radii of stars Q. Shortly 
afterwards it was recognized that a similar procedure can also be applied to pions 
emitted in high energy reactions to measure the size of the pion emitting source 
However, it is clear that the different conditions prevailing in high-energy or 
heavy-ion reactions ask for certain modifications of the original formalism. The 
table 1 lists the relevant parameters encountered in the study of stars or hot nu- 
clear matter. The correlation areas Af2 are given by the solid angle necessary to 



Table 1: Parameters relevant to stars and nuclei. 





star (photon) 


nucleus (pion) 


source size R 

momentum p 
correlation area Af2 = ixj (pR) 2 
de Broglie wavelength A = h/p 
coherence length A 


10 y m 
2 eV/c 
3 -lCT 32 sr 
6 -1CT 7 m 
CTfs3m 


10" 14 m 
2 -10 8 eV/c 
3 -lO" 2 sr 
6 -lO" 15 m 
A 2 /AA« 10" 14 m 



observe the speckle pattern in the emission of a large number of identical bosons as 
displayed in Fig.l of ref. [Sj. Their difference by 30. order of magnitudes is compen- 
sated by the equivalent difference in the distance between the star, respectively the 
nucleus and the interferometer. The large difference between the source size R and 
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the coherence length A in case of a star ensures that the star can be treated as a 
chaotic radiator. In general one finds for stars R S> A 3> A and the emitted photons 
may be described by plane waves. In case of the nucleus one finds R ~ A ~ A, and 
wavepackets are the more appropriate description of the emitted pions. 



2 Wavepacket Formalism 



The width of a Gaussian wavepacket at time t=0 is given by o"o which from the 
principles of quantum mechanics is related to the pion coherence length: do ~ 
A/47T ~ 1 fm.At later times t the free, i.e. non-interacting wavepacket in coordinate 
space is j|] 

¥(r,t) = (2^ 2 )~ ! exp jl (Pr - Et) - {v j , (1) 

where P, E are the centre momentum , energy of the wavepacket, R(i) is its spatial 
coordinate which changes with time according to R(t) = R + (P/m) • t, and the 
dispersion of the width is given by 



Similarly the wavepacket may be described in momentum space, where it turns out 
to be stationary: 



2 X (P-P) 2 



M^^rHJjafP-P)^)-^! (3) 



with dp = h/2ao. It is obvious that the wavepacket representation of the pions 
fulfills Heisenbergs uncertainty relation 

ft 

^ > g • ( 4 ) 

In the following we will restrict our discussion to SIS energies, i.e. the like- 
charge pion multiplicities are in central collisions of very heavy systems ~ 20, 
and multi pion effects can be neglected. We parametrize the pion source function 
by g(R,P) =p(R)f(P) with 

p(R) = ( 7T R*y l exp\-^\ (5) 
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Rl 

p2 | 



f(P) = (27rmTriexp{- — j ((5) 



and with the temperature T. The single pion observables in spatial and momen- 
tum space are then easy to calculate and yield the identical forms as in eqs.(5,6) 
except that the relevant parameters have to be replaced by R% = R^ + 2o"o and 
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Figure 1: Left panel: The exact and 1. order correlation functions. Right panel: The 
dependence of the chaoticity A on the localization cr . 



'-eff 



T + T qm where T qm 



n? 



Amo. 



This implies that for small pion multiplici- 



ties the localization do of the pions in the source does not change the single pion 
distributions, but the relevant parameters are changed: The effective source radius 
is modified by the pion localization, and the effective temperature is modified by a 
quantum contribution which takes into account the zero point energy of the pions. 

Similarly the two pion correlation function C2(pi,P2) = v^(p^)T\(pi) ma y Re- 
calculated. This calculation is more involved than in case of the single pion ob- 
servables, the exact result was published in As an approximate expression one 
obtains 2 

C 2 (pi,p 2 ) = l + Aezpj-^q 2 j , q = pi-p 2 (7) 



with an effective radius parameter R 2 



eff 



R s+ 2a 0^JJ 



The chaoticity parameter A is a function of the ratios R 2 s /{2aQ 
be expanded into the series 



so that R s < R e ff < R s - 
and T/T e fj and can 



A = 1 + 



oo 

E 

k=l 



(-)' 




2al 



1 - k- 



z 



eff, 



The Fig.l displays on the left side the exact and first order (A 



(8) 

1) results for 62(9) 
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using specific values for R s , <jq and T, and on the right side the dependence of A 
on do is shown for 3 different values of R s and the temperature T = 50 MeV. From 
this dependence we infer that A approaches only when ao becomes very large and 
T goes to zero. 

The wavepacket formalism for pions introduces a new parameter ao which was 
called localization because it determines the probability to localize the pion within 
a certain volume which is usually assumed to be within the hot nuclear matter. 
Naively the value of a$ is therefore experimentally bounded by the two limits R e ff = 
V2<jq and T e jf = T qm which yields 0.8 fm < do < 4 fm under SIS conditions. 
However it is more appropriate to relate ao to the coherence length A of the pion, 
which in the case that the pion is emitted in the decay of the A(1232) resonance 
is coupled to the lifetime of that resonance in hot nuclear matter. Alternatively 
one may also consider the mean free path of pions in hot nuclear matter as the 
appropriate quantity. We have used this latter conjecture in our calculations and 
have assumed ao = 1.8 fm. This implies a quantum contribution of T qm = 21 MeV 
to the effective temperature of pions. The plane wave limits are obtained in two 
ways: ao — ► oo (sharp momentum states) yields C2(pi,P2) = 1 + <5(pi — P2), in 
contradiction to experiment, ao — > (sharp position states) yields an infinite zero 
point energy, also in contradiction to experiment. This again demonstrates the 
inherent difficulties to describe the HBT interferometry with pions in a plane wave 
formalism. The wavepacket formalism was applied in other recent papers to the 
HBT interferometry with pions 0, [|, ||, @, §. In cases where ao was specified its 
value was set to ao ~ 0.7 fm which yields T qm ~ 100 MeV, a value which appears 
to be too high for the SIS energy range. 

3 The Coulomb Residual Interactions 

HBT interferometric studies in relativistic heavy-ion reactions are usually performed 
with charged pions. In this case the pions after emission cannot be described by 
free wavepackets but experience residual interactions from the charge of the source 
and from their own charges. As a result the centre motion of the wavepacket is 
modified and the wavepacket gets deformed. There are two ways to study these 
effects with respect to HBT interferometry, both ways are based on the time depen- 
dent Schrodinger equation: The multiconfigurational method |J, and the molecular 
dynamic method Q. 

3.1 Coulomb Interaction with the Source 

The second method was used to study the modifications of the two-pion correlation 
function when both pions interact with the source via the Coulomb force. Depending 
on the charge of the source the width of the correlation function is decreased for 
(7r - ,7r - ) but increased for (7r + ,7r + ) pairs, yielding a larger effective source radius 
from the former pairs, and a smaller radius from the latter pairs. The situation 
is depicted in the left panel of Fig. 2 for a source with radius R s = 5.5 fm and 
temperature T = 50 MeV. Surprisingly the increase respectively decrease are not 
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Z of source Zof 



Figure 2: Left panel: The dependence of the effective source radii from (tt~,tt~) and 
(7r + ,7r + ) pairs on the charge Z of the source. Right panel: The escape times t esc for 7r~ 
and 7r + as functions of Z. 



of same size, but stronger for the (tt~ ,tt~) pairs. The reason for this asymmetry is 
found in the right panel of Fig. 2 where the average time is shown which a charged 
pion spends within the escape radius R esc = 20 fm of the source. In case of the 
(tt + ,tt + ) pairs the correlation signal develops to a large fraction after the ir + have left 
the interaction region, and it is therefore only little influenced by this interaction. 
In case of the (ir~ ,tt~ ) pairs the times the tt~ spend inside the interaction region are 
longer and therefore this interaction disturbs the correlation signal more strongly. 
The asymmetry in the radius distortions is a direct proof that the appearance of 
the correlation does not occur instantaneously but is a process in time. 



3.2 Coulomb Interaction between Charged Pions 

The first and the second methods were used to study the effects of their mutual 
Coulomb interaction onto the two-pion correlation function. The results were sim- 
ilar: The Coulomb interaction between like-charge pions alters the correlation by 
such a small amount that it becomes unobservable under normal experimental con- 
ditions ||. This conclusion is in apparent contradiction to experimental results from 
(7r + ,7r - ) pairs [10], we will therefore focus our discussion on this type of pairs. 



The reason why the suppression of the correlation signal for very small values of q 
is not observed for like-charge pion pairs is connected to the zero point energy Eq = 
|T gm of the pion which is almost two orders of magnitude larger than the Coulomb 
energy between two pions with localization ctq = 1.8 fm. This is particularly easy 
to see in the case of two pions with different charges where this Coulomb energy 
attains its maximum value at complete overlap and is, for particles with charge 
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ure 3: Left panel: The time development of the relative wavepacket of unlike-charge 
with charge Z — 1. Right panel: The same for pairs with charge Z = 20. 



number Z, of size Ecoul 



Z 2 e 2 
\Ecoul 



One may define the parameter [11 1 



o 



^-Z mcro 



3vW 



(9) 



the value of which determines whether or not the Coulomb interaction becomes 
discernible. For a$ = 1.8 fm and Z = 1 one obtains £ = 1.4 • 10 -2 , but £ = 5.6 in 
case of a hypothetical charge number Z = 20. The time development of the pion 
wavepacket in relative coordinates is shown in Fig. 3 for both cases, the contour 
lines give the shape of the wavepacket, the dark curves the movement of its centre. 
It is evident that in the first case the wavepacket movement and its deformation 
is not modified by the Coulomb interaction whereas it is strongly modified in the 
second case with respect to the centre as well as to the deformation. In order 
to observe the Coulomb interaction the £ parameter has to be of order 1. This 
condition cannot only be met by an increase of the charge number Z but also by 
increasing the reduced mass of the particle pair or by increasing the localization uq. 
Finally one may also test the validity of the present procedure by approaching the 
classical limit h — > 0. These dependences were studied in detail, as examples the 
Fig. 4 displays in the left panel the (tt + ,7r - ) correlation functions for three different 
values of cjq and for the case fi — > ft/10. In this latter case the calculated result 
almost agrees with the classical expectation shown by the continuous curve. The 
dependence on o"o is as predicted, as a cross check, shown in the right panel, we 
have also performed the calculation for oq = 1.8 fm using the second method. 

We conclude that the Coulomb interaction only modifies the correlation function 
at small values of q when the £ parameter is of order one. Whereas this is not to 
be expected at SIS energies in case of pion pairs, it can occur for such pairs at 
much higher energies. At high energies pions are also emitted from the decay of 
longdiving resonances which were produced by the nucleus - nucleus collision. Such 
pions have a large coherence length with the related extended localization. The 
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Figure 4: The correlation function of (7r + ,7r~) pairs with different localizations &q. Also 
shown is the classical result (full curve) and its limit H — > h/10. The results in the left 
panel are from method 1, in the right panel from method 2 (see text). 



strength of the Coulomb distortion depends on the relative abundance of pions from 
the interaction zone to those from the long-living resonances, only the former will 
contribute to the observation of the HBT signal. The strength of the HBT signal 
at q — > also depends on the pion localization and on the source temperature, 
dependences which experimentally have not been explored up to now. Our results 
illustrate the important role the pion coherence length plays in describing the HBT 
interferometry with pions. 
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